New Relativistic Effects in the Dynamics of Nonlinear Hydrodynamical Waves 
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In Newtonian and relativistic hydrodynamics the Riemann problem consists of calculating the evolution of a 
fluid which is initially characterized by two states having different values of uniform rest-mass density, pressure 
and velocity. When the fluid is allowed to relax, one of three possible wave-patterns is produced, corresponding 
to the propagation in opposite directions of two nonlinear hydrodynamical waves. New effects emerge in a 
special relativistic Riemann problem when velocities tangential to the initial discontinuity surface are present. 
We show that a smooth transition from one wave-pattern to another can be produced by varying the initial 
tangential velocities while otherwise maintaining the initial states unmodified. These special relativistic effects 
are produced by the coupling through the relativistic Lorentz factors and do not have a Newtonian counterpart. 

PACS Numbers: 47.40.N, 47.35, 47.85.D 



More than a hundred years ago, Riemann addressed the 
problem of the one-dimensional time evolution of a perfect 
fluid which, at some given initial time t = 0, is composed of 
two adjacent states characterized by different values of uni- 
form velocity, pressure and density [[[]]. Since then this prob- 
lem has been referred to as the "Riemann problem" and is the 
prototype of the initial value problem for hyperbolic systems 
of partial differential equations with discontinuous initial con- 
ditions. The conclusion reached by Riemann in Newtonian 
hydrodynamics is that the one-dimensional flow that develops 
when the barrier separating the initial "left" (or 1) and "right" 
(or 2) states is removed, will allow for four different and dis- 
tinct solutions. All of the solutions are composed of nonlinear 
waves, in the form of either shock waves or rarefaction waves, 
that propagate in opposite directions and join the two unper- 
turbed left and right states. As a result, three different "wave- 
patterns" can be produced corresponding to two shock waves 
(2S), a shock and rarefaction wave (SR) and two rarefaction 
waves (2i?), respectively. Schematically, the fluid solution for 
any t > can therefore be represented as 

1W^3C3'W_2, (1) 

where W denotes a shock or a rarefaction wave that prop- 
agates towards the left (<— ) or towards the right (— ►) of 
the initial discontinuity, 1 and 2 are the known initial states 
with "state-vectors" consisting of fixed and independent val- 
ues of the pressure, rest-mass density and velocity U1.2 = 
(p,p,f)i,2. Similarly, 3 and 3' represent the new hydrody- 
namic states that form behind the two waves propagating in 
opposite directions. The two nonlinear waves W<- and 
are separated by a contact discontinuity C across which the 
rest-mass density is discontinuous, but not the velocity nor 
the pressure, i.e. p^ ^ py, p^ = py, v$ = V31. These three 
quantities fully describe the states of the fluid in the flow re- 
gions 3 and 3' behind the nonlinear waves. The formal so- 
lution of the Riemann problem consists then in determining 
the velocity, pressure and rest-mass density in the new states 



3 and 3' as well as calculating, at any time t > 0, the positions 
of the waves separating the four states. 

The solution of the Riemann problem attracted a wider in- 
terest when it was realized that its numerical solution could 
be implemented as building blocks in the construction of nu- 
merical Godunov-type methods for the accurate solution of 
the hydrodynamical equations [Q], In such methods, the com- 
putational domain is discretized and each interface between 
two adjacent grid-zones is used to construct the initial left and 
right states of a "local" Riemann problem. The evolution of 
the hydrodynamical equations is then obtained through the so- 
lution of the sequence of local Riemann problems set up at the 
interfaces between successive grid-zones and basically con- 
sists of determining the fluid pressure in regions 3 and 3' (see 
ref. [|3|] for a recent review). 

The extension of the one-dimensional Riemann problem 
to relativistic fluid velocities has been performed by [|J for 
the "shock-tube" problem and by [|5|] for the general problem. 
Besides the obvious additional complications introduced by 
special relativity, the one-dimensional Riemann problem in 
relativistic hydrodynamics does not show qualitative differ- 
ences from its Newtonian counterpart. In a recent paper [Q], 
we have introduced a new approach to the solution of an ex- 
act relativistic one-dimensional Riemann problem which fo- 
cuses on the relativistic invariant expression for the relative 
velocity between the two unperturbed initial states (vf 2 )o = 
[(vf — uf)/(l — w i u f)]o> where we have assumed that the 
discontinuity is initially located on an x = const, surface and 
where v x represents the component of the 4-velocity normal 
to this surface. This new approach introduces great simplifica- 
tions in the logical formulation of the Riemann problem and 
has been essential in revealing these new relativistic effects. 
For this reason we will briefly review it in what follows. 

By construction, (vf 2 )o measures the relativistic jump of 
the fluid velocity normal to the discontinuity surface. Within 
this new approach, the solution to the relativistic Riemann 
problem is found after the pressure in the region between the 
two nonlinear waves has been calculated as the root of the 
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nonlinear equation [g] 



K 2 )o - o . 



(2) 



The notation used in equation (|2j) needs an additional com- 
ment. It is should be emphasized, in fact, that p 3 itself de- 
pends on the thermodynamic al properties of the initial states 



so that vf 2 
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(pi , Pi, P2i Pi) • Note also that vf 2 (p 3 ) 



does 



not depend on the value of the initial tangential velocities and 
has a functional form that is different for each of the three 
possible wave-patterns that might result from the decay of the 
initial discontinuity. In ref. [^] it was then shown that the 
wave-pattern produced in a one-dimensional relativistic Rie- 
mann problem by the decay of the discontinuity can be pre- 
dicted from knowledge of the initial data. The possibility of 
predicting the wave-pattern from the initial conditions applies 
also to Newtonian hydrodynamics [Q] and represents an im- 
portant advantage since it allows one to deduce in advance 
which set of equations to use for the solution of the Riemann 
problem. From a mathematical point of view, the validity of 
the approach discussed in ref. [||] is based on the proof that 
the three branches of the function vf 2 — vf 2 (p 3 ), correspond- 
ing to the three possible wave-patterns, are monotonically in- 
creasing with p 3 . Furthermore, it was shown that the three 
different branches always join smoothly through specific val- 
ues of vf 2 (ps) denoted respectively as (vf 2 ) 2S and (vf 2 )sR- 
(In ref. |Q] another limiting value, (vf 2 ) 2R , was found but this 
is not relevant in the present discussion and will not be con- 
sidered here.) 

Exploiting these properties, an improved exact Riemann 
solver can be built in which the comparison of (vf 2 )o with 
the relevant limiting values {v'f 2 ) 2S and (vf 2 ) SR constructed 
from the initial conditions allows to determine, prior to solv- 
ing equation (^), both the wave-pattern produced and the func- 
tional form of vf 2 (ps) to be used in equation (Q) [^]. In prac- 
tice, this approach leads to a straightforward numerical im- 
plementation and a more efficient numerical algorithm with 
computational costs reduced by up to 30%. 

In generic multidimensional flows, the Riemann problem 
can still be cast in terms of a planar discontinuity surface 
across which a normal flow with velocity v x takes place and 
is responsible for the mass transfer. In general, however, the 
flow will also have components tangential to the discontinuity 
surface, which we will refer to as v* = \{v v ) 2 + {v z ) 2 ] 1 ^ 2 . The 
latter can also be discontinuous across C possibly leading to a 
Kelvin-Helmoltz instability there. In relativistic regimes, the 
tangential component cannot be removed by a single Lorentz 
boost and is constrained to satisfy \v l \ < [1 — (v x ) 2 } 1 ^ 2 . 

A first study of a multidimensional relativistic Riemann 
problem involving two strong rarefactions was carried out in 
ref. [ 3 ] and has been extended to a more general treatment in 
ref. [9] where it was pointed out that the tangential velocities 
in multidimensional flows introduce differences between the 
Newtonian and the relativistic Riemann problems. In New- 
tonian hydrodynamics, in fact, the solution to the Riemann 
problem does not depend on the tangential component of the 



flow. In relativistic hydrodynamics, on the other hand, the 
different regions of the flow are coupled, both in the veloci- 
ties and in the specific enthalpy, through the Lorentz factors 
Wi = [1 - (uf) 2 - (v*) 2 ]~ 1/2 , where i = 1,2. This is an 
important difference whose consequences are much harder to 
assess. 
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FIG. 1. Relative velocity as a function of the pressure in region 
3. The continuous line refers to the case of zero tangential velocities 
while the dashed and dotted ones refer to cases in which the tangen- 
tial velocities are nonzero (We use units in which c = 1.). The filled 
circles represent the joins between the solutions consisting of two 
shock waves (2S), of one shock and one rarefaction wave (SR) and 
of two rarefaction waves (2R). 

To discuss the new effects that emerge in a multidimen- 
sional and relativistic Riemann problem we have plotted in 
Fig. [j] the function vf 2 for different values of the pressure at 
the contact discontinuity p 3 . The curves shown refer to the 
case of zero tangential velocities (continuous line) and to the 
cases where the tangential velocities are nonzero (dashed and 
dotted lines), respectively (Note that the curves do not depend 
on the sign of v*.). In both cases the relative velocities are 
calculated for a perfect fluid with polytropic equation of state 
p = kp' f , with a standard set of the initial values of pressure 
and rest-mass density [10], i.e. p\ = 1.0, P2 = 0.1, p\ — 
1.0, p2 — 0.125 (we have here chosen 7 = 5/3). The filled 
dots in each curve indicate the locations where the branches 
representing the three wave-patterns (2S, SR, and 2R) merge. 
The values of the relative normal-velocities where this hap- 
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pens are just the values of 
above. 

A rapid examination of Fig. [j] is sufficient to appreciate that 
the presence of tangential velocities can introduce new qual- 
itative differences in a relativistic Riemann problem. When 
tangential velocities are present, in fact, the relative normal- 
velocity vf 2 is a function of p 3 but, through the coupling in- 
troduced by the Lorentz factors, also of v\ and v\. This can 
introduce a new relativistic effect and shift the solution from 
one wave-pattern to another. To better appreciate this, let us 
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restrict attention to a simpler situation in which only one of 
the two initial tangential velocities is allowed to vary. 
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FIG. 2. Solution of the Riemann problem with initial conditions: 
Pi — 1.0, P2 = 0.1, pi — 1.0, p2 = 0.125. The initial velocities 
are v\ = 0.0 = v\ and vf = 0.5, «f = 0.0. 

Suppose now that the normal velocities are chosen to be 
vf = 0.5, v 2 = 0.0, and that there are no tangential ve- 
locities. In this case, (vf 2 )o = 0.5 and Fig. [I] shows that the 
solution to the Riemann problem falls in the SR branch, hence 
producing a wave-pattern consisting of a shock and a rarefac- 
tion wave moving in opposite directions. This is presented in 
Fig. ^ which shows the solution of the Riemann problem at 
a time t > for the pressure, the rest-mass density and the 
velocity. However, if we now maintain the same initial con- 
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FIG. 3. As in Fig. | but with v\ = 0.0, v\ = 0.9, vf = 0.5, 
uf = 0.0. 

ditions but allow for nonzero tangential velocities in state 2, 
Fig. [I] also shows that the solution to the Riemann problem 
can fall in the 25 branch (cf. dotted line), hence producing 
a wave-pattern consisting of two shock waves moving in op- 
posite directions. This is presented in Fig. [5] which shows the 
solution of the same Riemann problem as in Fig. but with 



the initial tangential velocities being v\ = 0.0 and v 2 = 0.9. 
Note that except for the tangential velocities, the solutions in 
Figs. U and || have the same initial state-vectors but different 
intermediate ones (i.e. p%, p%, and i>f). 

The Riemann problems shown in Figs. || and || are only two 
possible examples but they are useful ones since they show 
that through the coupling among the different states intro- 
duced by the Lorentz factors, a difference in the tangential 
velocities can produce a smooth transition from one wave- 
pattern to another while maintaining the initial states unmod- 
ified. Interestingly, this transition does not need to always 
produce a solution consisting of two shock waves. 
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FIG. 4. As in Fig. | but with v\ = 0.0 = v\, vf = 0.0, uf = 0.5. 

Suppose, in fact, that the normal velocities are now chosen to 
be vf = 0.0, t>2 = 0.5. We can then repeat the considera- 
tions made above and start examining the wave-pattern pro- 
duced when there are zero tangential velocities. In this new 
setup, (1*12)0 = —0.5 and Fig. [j] shows that the solution to 
the Riemann problem still falls in the SR branch (cf. dashed 
line), with the corresponding solution at a time t > being 
presented in Fig. f|. Note that the wave-patterns in Fig. ||and^ 
both consist of a shock and a rarefaction wave, but have alter- 
nate initial normal velocities. 

When nonzero tangential velocities are now considered in 
state 1, Fig.|l] shows that {vf 2 )o can fall in the 2R branch, 
hence producing a wave-pattern consisting of two rarefaction 
waves moving in opposite directions. This is again shown 
in Fig. U where we have chosen initial tangential velocities 
v\ = 0.9 and v\ = 0.0. In this case too and except for the 
tangential velocities, the solutions in Figs. |] and |] have the 
same initial state-vectors but different intermediate ones. 

As mentioned above, the appearance of these new relativis- 
tic effects is related to the behaviour of the function vf 2 for 
different values of the initial tangential velocities and in par- 
ticular to how the three branches composing the curve change 
under variation of v\ 2 . As a result, the occurrence of these 
effects can be recast into the study of the functional behaviour 
of the values of vf 2 delimiting the different branches when the 
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FIG. 5. As in Fig. | but with v\ = 0.9, v\ = 0.0, vf = 0.0, 
uf = 0.5. 

tangential velocities are varied within the allowed ranges. A 
detailed mathematical treatment of this problem will be pre- 
sented in a separate paper [|ll]] but it is sufficient to point out in 
what follows the main results. Firstly, the dependence of the 
limiting values for the SR and 25 branches on the tangential 
velocities can be expressed as (cf. Fig. [j]) 

(S?a) M = (0? 2 ) M («i) . («?a) as = (5?a) a >£) ■ (3) 

Secondly, these values converge to zero when the tangential 
velocities are allowed to reach their maximum values, i.e. 

Jim (vf 2 ) SR = = Jim (vf 2 ) 2S . (4) 

VVi^oo W2— >oc 

Expressions (||) indicate that for tangential velocities assum- 
ing increasingly larger values, the SR branch of the vf 2 curve 
spans a progressively smaller interval of relative normal- 
velocities. When the tangential velocities reach their asymp- 
totic values, the SR branch reduces to a point. As a result, a 
wave-pattern consisting of a shock and a rarefaction wave is 
generically disfavoured by increasing tangential velocities. 

The effects discussed so far have a purely special relativis- 
tic origin and might conflict with our physical intuition. How- 
ever, a behaviour similar to the one reported here is found 
in the well-known relativistic transverse-Doppler effect, in 
which the wavelength of a photon received from a source 
moving at relativistic speeds changes also if the source has 
a velocity component orthogonal to the direction of emission 
of the photon [Q. In this case too, a Lorentz factor including 
the transverse velocity is responsible for the effect. 

A couple of comments are worth mention. Firstly, there ex- 
ists a set of initial conditions for which these new relativistic 
effects will not occur. These initial conditions are those of the 
classic "shock-tube" problem, in which vf = = v% . In this 
case, in fact, {vf 2 )o = and, because of the limits (Q), the so- 
lution of the Riemann problem will be given by a wave-pattern 
consisting of a shock and a rarefaction wave, independently of 



the values of the tangential velocities. Secondly, while the sta- 
bility of shock fronts has been studied in the past (see [[l3|] and 
references therein), little is still know about the stability prop- 
erties of shock fronts across which the tangential velocities are 
discontinuous. It is plausible that a nonzero mass flux makes 
these fronts stable with the respect to the Kelvin-Helmoltz in- 
stability, but more work needs to be done in the investigation 
of the stability of generic shocks fronts. 

The results reported in this letter can be used in astrophysi- 
cal scenarios, such as those involving relativistic jets or 7— ray 
bursts, in which nonlinear hydrodynamical waves with large 
Lorentz factors and complex multidimensional flows are ex- 
pected [53]. Furthermore, both the approach to the solution 
of the Riemann problem and the new relativistic effects dis- 
cussed here could be relevant in the implementation of spe- 
cial relativistic exact Riemann solvers in multidimensional 
codes solving the hydrodynamical equations in curved space- 
times [|,[l5[]. 
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